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ABSTRACT 

We report the discovery of a confirmed supernova (SN) and a supernova-candidate in near-infrared 
images from the ALTAIR/NIRI adaptive optics system on the Gemini-North Telescope and NICMOS 
on the Hubble Space Telescope. The Gemini images were obtained as part of a near-infrared K- 
band search for highly-obscured SNe in the nuclear regions of luminous infrared galaxies. SN 2008cs 
apparent in the Gemini images is the first SN discovered using laser guide star adaptive optics. It 
is located at 1500 pc projected distance from the nucleus of the luminous infrared galaxy IRAS 
17138-1017. The SN luminosity, JHK colors and light curve are consistent with a core-collapse event 
suffering from a very high host galaxy extinction of 15.7 ± 0.8 magnitudes in U-band which is to our 
knowledge the highest yet measured for a SN. The core-collapse nature of SN 2008cs is confirmed 
by its radio detection at 22.4 GHz using our Very Large Array observations 28 days after the SN 
discovery, indicating a prominent interaction of the SN ejecta with the circumstellar medium. An 
unconfirmed SN apparent in the NICMOS images from 2004 is located in the same galaxy at 660 pc 
projected distance from the nucleus and has a lower extinction. 

Subject headings: supernovae: individual(SN 2008cs) - galaxies: starburst - galaxies:individual(IRAS 
17138-1017) - infrared: galaxies - instrumentation: adaptive optics 



1. INTRODUCTION 

A large fraction of the massive star formation at high-z, 
e.g., 70% at z = 1 (Le Floc'h et al. 2005), took place in 
luminous (Lir, > 10 11 L©) and ultraluminous (Lir, > 10 12 
L ) infrared (IR) galaxies (LIRGs and ULIRGs, respec- 
tively), and their high star formation rates (SFRs) can 
be expected to result in core-collapse supernova (CCSN) 
rates a couple of orders of magnitude higher than in ordi- 
nary field galaxies. Most of the SNe occurring in LIRGs 
and ULIRGs are likely to be obscured by large amounts 
of dust in the nuclear starburst environment; however, in 
the near-IR the extinction is strongly reduced. Because 
of the likely concentration of most of the SNe within the 
innermost nuclear regions a high spatial resolution is cru- 
cial for their detection. This can be achieved with space- 
based imaging (e.g., Cresci et al. 2007, Colina et al. 
2007), imaging from Antarctica (see Burton et al. 2005) 
or ground-based adaptive optics (AO) imaging provid- 
ing resolutions 5 times better than under typical natural 
seeing conditions. The potential of the current 8m-class 
telescopes equipped with AO was recently demonstrated 
by the discovery of SN 2004ip in the LIRG nuclear re- 
gions using the NACO AO system with a natural guide 
star on the Very Large Telescope (Mattila et al. 2007). 

In this letter we report a discovery of two SNe in 
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IRAS 17138-1017 which is a LIRG (L IR = 2.4 x 10 11 
L Q , Sanders et al. 2003) at a redshift of z = 0.0173 
(Shier & Fischer 1998) corresponding to a distance of 
75 Mpc (H = 70 km s^Mpc" 1 , ft A =0.7 and fi M =0.3). 
The SNe were detected as a result of our X-band search 
for highly-obscured CCSNe in a sample of nearby LIRGs 
using the Gemini-North Telescope. Even though IRAS 
17138-1017 has a high expected SN rate, only one SN 
(the Type la SN 2002bw, with a likely low extinction; Li 
2002, Matheson et al. 2002) has been detected in this 
galaxy before. 

2. OBSERVATIONS AND RESULTS 

2.1. IR observations 

IRAS 17138-1017 was observed with iKT/NICMOS 
as part of a program to observe a complete sample of 
LIRGs (Alonso-Herrero et al. 2006) with the NIC2 cam- 
era (0.076" pixel" 1 , FWHM - 0.1-0.2") on 2004 Septem- 
ber 23.0 UT (see Table 1). Observations were conducted 
with the F110W and F160W broadband filters and with 
the F187N narrowband filter. Details of the data reduc- 
tion and photometric calibration of the NICMOS images 
can be found in Alonso-Herrero et al. (2006) . The LIRG 
was next observed with the Gemini-North Telescope 7 in 
the if -band using the Near-InfraRed Imager (NIRI) with 
the ALTAIR laser guide star AO system (0.022" pixel" 1 , 
FWHM - 0.1") on 2008 April 21.6 UT (program GN- 
2008A-Q-38, PI: S. Ryder). The NIRI images were re- 
duced using the NIRI package in iraf. The target frames 

7 Based on observations obtained at the Gemini Observatory, 
which is operated by the Association of Universities for Research 
in Astronomy, Inc., under a cooperative agreement with the NSF 
on behalf of the Gemini partnership: the National Science Founda- 
tion (USA), the Science and Technology Facilities Council (UK), 
the National Research Council (Canada), CONIC YT (Chile), the 
Australian Research Council (Australia), Ministrio da Cincia e Tec- 
nologia (Brazil) and SECYT (Argentina). 
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were de-dithered and median combined in iraf using a 
bright field star as a reference point. Comparison be- 
tween the HST and Gemini observations revealed two 
new point sources, one at each epoch. The object ap- 
parent only in the HST images is located at 1.9" (or 
660 pc) and the one seen only in the Gemini images at 
4.2" (or 1500 pc) projected distance from the near-IR nu- 
cleus of IRAS 17138-1017. Both objects were reported to 
the Central Bureau for Astronomical Telegrams (CBAT) 
immediately after discovery and the more recent one was 
later assigned the designation SN 2008cs (Kankare et al. 
2008, Perez- Torres et al. 2008). Follow-up observations 
of SN 2008cs were obtained on May 20.5 UT with NIRI 
in the J, H, and K bands without AO (0.117" pixel -1 , 
FWHM - 0.5"), and with ALTAIR/NIRI on May 27.5, 
June 25.5, and September 20.2 - 21.2 UT (see Fig. 1) us- 
ing the AO. Follow-up observations were also obtained 
using the Nordic Optical Telescope (NOT) with the 
NOTCam instrument (0.234" pixel" 1 , FWHM - 0.7") in 
the J, H, and K s bands on June 18.1 UT (see Table 1). 

We used 35 stars from the 2MASS catalog to de- 
rive an accurate World Coordinate System (WCS) for 
an HST/ACS /-band image of IRAS 17138-1017 ob- 
tained from the HST Science Archive. The Gemini/NIRI 
if-band (AO) image from June 25 was aligned to the 
ACS image using seven pointlike sources common be- 
tween the frames and a simple geometric transforma- 
tion including shifts in x and y, and the same scale 
factor and rotation for both x and y. This yielded 
17 h 16 m 35.855 s and -10°20'42.98" for SN 2008cs with 
0.09" and 0.07" uncertainty in R.A. and Deck respec- 
tively. The ffST/NICMOS F160W image was aligned 
to the ACS image using 6 point sources and this yielded 
R.A. = 17 h 16 m 35.900 s and Deck = -10°20'37.85" for 
the 2004 event. 

2.2. SN 2008cs photometry and extinction 

The photometry of SN 2008cs was measured from the 
images using the SNOOPY 8 package in iraf. Suitable 
field stars were selected to fit a point-spread function 
(PSF) for each image. A polynomial surface was fitted 
to a background region around the object position (ex- 
cluding the innermost region around the object) and af- 
ter removing the local background the PSF was fitted to 
the object and a residual image was produced after sub- 
tracting the PSF. Error estimates for the measurements 
were obtained by simulating and PSF fitting nine artifi- 
cial point sources around the SN in the residual images. 
For the non-AO Gemini/NIRI and NOT images the pho- 
tometry was calibrated using six 2MASS stars within the 
image field of view (FOV). The AO Gemini/NIRI images 
with small FOV were photometrically calibrated to the 
non-AO images using PSF photometry of three field stars 
available in both the AO and non-AO images (see Fig. 1). 
Assuming that these field stars have a similar AO PSF 
as the SN, this procedure also corrects for the flux lost in 
the PSF halo around the diffraction limited core in the 
AO images. Finally, all the magnitudes were aperture 
corrected using three bright isolated stars in the non-AO 
images to an aperture of 3" in radius. For the complete 

8 SNOOPY, originally presented in Patat (1996), has been im- 
plemented in IRAF by E. Cappellaro. The package is based on 
DAOPHOT, but optimised for SN magnitude measurements. 



TABLE 1 

Log of the observations and supernova photometry. 
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Sept. 23.0 


INlGMOb, r 11UW 


3 X 


80 s 


18.28 ± 0.U6 


Sept. 23.0 


JNlOMUb, r IbUW 


3 x 


80 s 


17. UU ± U.U5 


Sept. 23.1 


IMLiMUo, rlo/IN 


3 x 


288 s 


lb. do ± U.Ub 


2008 








SN 2008cs 


April 21.6 


A 1 PITTTn y\TTTST 7S 

ALTAIR/NIRI, K 


8 x 


30 s 


16.65 ± 0.06 


May 20.5 


NIRI, K 


9 x 


30 s 


16.22 ± 0.37 


May 20.5 


NIRI, H 


7 x 30 s 


17.40 ± 0.46 


May 20.5 


NIRI, J 


9 x 


30 s 


19.43 ± 0.50 


ivlay z ( .o 


ATTATP /T\n"DT If 
ALjlAlrx/lNlixl, A 


8 x 


30 s 


10. it ± U.UD 


May 27.5 


ALTAIR/NIRI, H 


8 x 


30 s 


16.97 ±0.06 


June 18.1 


NOTCam, K s 


18 x 


60 s 


15.99 ±0.61 


June 18.1 


NOTCam, H 


18 x 


60 s 


16.92 ± 0.35 


June 18.1 


NOTCam, J 


18 x 


60 s 


19.58 ±0.65 


June 25.5 


ALTAIR/NIRI, K 


8 x 


30 s 


15.82 ± 0.10 


June 25.5 


ALTAIR/NIRI, H 


8 x 


30 s 


17.17 ±0.05 


June 25.5 


ALTAIR/NIRI, J 


8 x 


30 s 


19.02 ± 0.14 


Sept. 20.2 


ALTAIR/NIRI, K 


8 x 


30 s 


17.30 ± 0.09 


Sept. 20.2 


ALTAIR/NIRI, H 


8 x 


30 s 


19.11 ±0.09 


Sept. 21.2 


ALTAIR/NIRI, J 


8 x 


30 s 


> 21.4 



list of apparent SN magnitudes see Table 1. 

Assuming a distance of 75 Mpc and a Galactic extinc- 
tion of Ay = 2.29 (Schlegel et al. 1998, but see also Arce 
& Goodman 1999) the absolute magnitude of the SN at 
the time of the second AO observation (when the SN is 
closest to the peak) becomes Mh = —17.81 and Mr = 
— 18.87 (before correcting for the host galaxy extinction). 
We compared these with the absolute peak magnitudes of 
typical Type II-P SNe 1999em and 2003hn (M H « -17.8 
and Mk ~ —18.1; Krisciunas et al. 2008) and the "or- 
dinary" (Mh = -18.58 and M K = -18.62) and "slowly- 
declining" (Mh = -19.95 and M K = -20.02) CCSN tem- 
plate light curves from Mattila & Meikle (2001). The ob- 
served SN color and peak absolute magnitude can only 
be accounted for by the "slowly-declining" type, and 
even then only with a substantial host galaxy extinction. 
Adopting these JHK templates and the extinction law 
of Rieke & Lebofsky (1985) for which A K = 0.112 x A v 
we obtain a host galaxy extinction towards SN 2008cs of 
1.67 — 1.85 mag in K. This corresponds to 14.9 — 16.5 mag 
in V (Ay = 17.2 — 18.8 total line of sight extinction) with 
the intrinsic absolute peak magnitude Mk = —20.9 ±0.1 
(statistical error only) occurring around 49 — 55 days af- 
ter the discovery (see Fig. 2 for an example fit). We note 
that in NICMOS images there is no star forming region 
detected at the position of SN 2008cs as expected due to 
the high extinction towards the SN (see Diaz-Santos et 
al. 2008). 

2.3. The CCSN rate of IRAS 17138-1017 

To estimate the SFR for IRAS 17138-1017 we compiled 
its 12-100 fim IRA S photometry and fitted the resulting 
spectral energy distribution (SED; see Fig. 3) with the 
starburst models of Efstathiou et al. (2000). The star- 
burst model provides a very good fit to the near- to far- 
IR SED and is consistent with a negligible contribution 
from an active galactic nucleus (AGN; see also DcPoy 
et al. 1988). The model assumes an exponentially de- 
caying starburst with a time constant of 30 Myr and an 
age of 57 Myr. The initial Ay of the molecular clouds 
that constitute the starburst is 100 mag. The SFR at 
the peak is 169 M Q yr _1 , and 76 M Q yr _1 when aver- 
aged over the duration of the starburst. We estimate a 
bolometric luminosity of 2.3 x 10 11 L for IRAS 17138- 
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Fig. 2. — An example fit for the SN 2008cs light curves with 
Ay(host) = 15.8. The SN is 0.9 mag brighter than the "slowly- 
declining" CCSN template light curves and the 7-f-band peak oc- 
cured at 53 days after the discovery, see text for details. 
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Fig. 1.— 14" X 14" subsections of IRAS 17138-1017 images taken with the HST /NICMOS F160W filter on September 23, 2004, and with 
Gemini ALTAIR/NIRI AO in #-band and J-band on June 25, 2008. The "historical" SN from 2004 and SN 2008cs are indicated in the 
NICMOS and Gemini images, respectively. The effects of high extinction in the SN 2008cs line-of-sight can also be seen by comparing 
the brightness of the SN to the field stars in the H- and J-bands. The three field stars indicated in the right-hand panel were used for 
photometric calibration. North is up and East is to the left in all the images. 

2.4. VLA observations 

We observed IRAS 17138-1017 on 2008 May 19.4 UT 
(program AP558, PI: M. Perez- Torres) at 8.4 and 22.4 
GHz, using the Very Large Array (VLA) in the C config- 
uration, aimed at detecting any radio emission from SN 
2008cs with rapid response exploratory VLA time. The 
observations lasted for 2.5 hr, with effective on-source 
times on the nucleus of IRAS 17138-1017 of -60 and -49 
minutes at 8.4 and 22.4 GHz, respectively. We used the 
source 1713-130 for phase and amplitude calibration pur- 
poses, and the quasar 3C 286 (1331+305) to set the ab- 
solute flux density scale. We then edited, calibrated, and 
imaged our 8.4 and 22.4 GHz VLA data, following stan- 
dard data reduction techniques implemented within the 
NRAO Astronomical Image Processing System (AIPS). 
Since the VLA was in C-configuration, the angular reso- 
lution of the 8.4 GHz observations did not allow the de- 
tection of any compact radio emission, apart from that of 
the nuclear region. But our 22.4 GHz (beam size ~ 0.9") 
observations clearly showed a local maximum of radio 
emission (£?„ = 440 ± 75/xJy, Perez- Torres et al. 2008) 
at R.A. = 17 h 16 m 35.853 s , Decl. = -10°20'43.05". The 
radio source appears coincident within the uncertainties 
(0.2" in each coordinate) with the position of the SN at 
near-IR wavelengths (offset of 0.03" and 0.07" in R.A. 
and Decl., respectively), which confirms that both the 
near-IR and radio emission come from the same region 
(see Fig 4). Future higher angular resolution observa- 
tions should also detect SN 2008cs at longer wavelengths. 

Type la SNe are not expected to be strong radio emit- 
ters and have not been detected at radio wavelengths 
(e.g., Panagia et al. 2006). Our detection of SN 2008cs 
at 22.4 GHz indicates a prominent interaction of the SN 
ejecta with the circumstellar medium, confirming that 
the SN was a core-collapse event. SN 2008cs is only the 
third SN after SNe 2000ft (Colina et al. 2001) and 2004ip 
(Perez- Torres et al. 2007) detected at both optical/IR 
and radio wavelengths in a LIRG/ULIRG host. 

2.5. The "historical" supernova from 2004 

We carried out the photometry of the 2004 event with 
circular apertures and subtracted the local background. 
The aperture corrections were done using synthetic PSFs 
generated with TinyTim (Krist et al. 1998). An as- 
teroid origin for the 2004 source was immediately ruled 
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Fig. 3.— The IR SED of IRAS 17138-1017 fitted with the 
starburst models of Efstathiou et al. (2000). No AGN contribution 
is needed to explain the SED. 



1017 which is close to its IR luminosity as expected for a 
LIRG. Adopting the average SFR and assuming CCSN 
progenitors between 8 and 50 Mq and a Salpeter initial 
mass function with cut-offs at 0.1 and 125M0, we obtain 
a CCSN rate of -0.53 yr^ 1 for IRAS 17138-1017. Using 
76 Mq yr -1 as an upper limit for the SFR averaged over 
0.5 Gyr we obtain an upper limit for its Type la SN rate 
(Sullivan et al. 2006) of -0.03 yr" 1 which is -20 times 
smaller than the estimated CCSN rate. 
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17 h 16 m 36.0 s 35.8 s 35.6 s 

Fig. 4.— if -band Gemini ALTAIR/NIRI image from June 25, 
2008 with 22.4 GHz VLA-C contours from May 19, 2008 overlaid. 

out since it is stationary compared to the field stars be- 
tween the first and the last target frame separated by 
~40 minutes. Also no known minor planets were near 
the source within a radius of 10' at the time of the obser- 
vation. A variable foreground star also provides a very 
improbable alternative explanation given the small FOV 
(19.8" x 19.8") of the HST /NICMOS frames. For exam- 
ple, the likelihood of catching a foreground nova would 
be extremely small given the Galactic nova rate of ~35 
yr- 1 (Shafter 1997). Based on this and the high CCSN 
rate of IRAS 17138-1017, we conclude that the 2004 tran- 
sient source was very likely a CCSN. From its (J-H) color 
and assuming a negligible difference between the F160W 
and the ground-based H filter and a conversion from the 
F110W to the ground-based J-band photometry (Origlia 
& Leitherer 2000) as well as the template CCSN light 
curves from Mattila & Meikle (2001) we can estimate a 
host galaxy extinction for the SN. We obtain Ay ~ 4 



mag if it was an "ordinary" CCSN near the peak and 
Ay ~ mag if it was a "slowly-declining" CCSN at ^4 
months from the peak. 

3. CONCLUSIONS 

Near-IR-discovered SNe in LIRGs and ULIRGs are not 
only important for SN rate estimates. They can also be 
used to estimate the amount and distribution of extinc- 
tion within their host galaxies, and follow-up observa- 
tions can provide us with a better understanding of the 
behavior of SNe within their dusty, high-density envi- 
ronments. To our knowledge the host galaxy extinction 
of A v = 15.7 ± 0.8 mag measured for SN 2008cs is the 
highest observed for a SN so far. Previous events with 
high extinctions include SN 200 ldb (Ay rj 5.5; Maiolino 
et al. 2002), SN 2002cv (Ay w 9; Di Paola et al. 2002, 
Elias-Rosa et al. 2008), SN 2002hh (Ay « 5; Pozzo et al. 
2006), SN 2004am (Ay w 5; Mattila et al., in prepara- 
tion), and SN 2004ip (Ay « 5 - 40; Mattila et al. 2007). 
Such high extinction events can have an important im- 
pact for the SN statistics when estimating the complete 
CCSN rates including also the optically obscured SNe. 
This will be essential when using CCSNe as probes of 
the SFR at both low- and high-z (e.g., Cappellaro et al. 
1999, Dahlen et al. 2004, Botticella et al. 2008, Melinder 
et al. 2008) with the aim of providing a new independent 
measurement of the cosmic star formation history. 
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